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We have measured the ^'"^Sm(n,a) cross section from 3 eV to 500 keV and performed an TZ- 
matrix analysis in the resolved region 700 eV) to extract a widths for 104 resonances. We 

computed strength functions from these resonance parameters and compared them to transmission 
coefficients calculated using optical model potentials similar to those employed as inputs to statistical 
model calculations. The statistical model often is used to predict cross sections and astrophysical 
reaction rates. Comparing resonance parameters rather than cross sections allows more direct tests 
of potentials used in the model and hence should offer greater insight into possible improvements. 
In particular, an improved a+nucleus potential is needed for applications in nuclear astrophysics. 
In addition to providing a more direct test of the a+nucleus potential, the a-width distributions 
show indications of non-statistical effects. 

PACS numbers: 



I. INTRODUCTION 



Reactions involving a particles and intermediate-to- 
heavy mass nuclides often can play an important role in 
nucleosynthesis occurring in massive stars at high tem- 
peratures and in explosive environments such as super- 
novae [1, 2]. For example, accurate rates for many (7,0;) 
reactions are needed for a better understanding of the 
nucleosynthesis of the neutron deficient A > 90 nuclides 
during the so-called p process. As important as these 
reactions are, there is scant experimental information 
on their rates because the cross sections are extremely 
small and many of the required "target" isotopes have 
very small natural abundances (and hence are very ex- 
pensive) or are radioactive; thus, direct measurements 
are very difficult. However, most of the required rates 
should be calculable to sufficient accuracy using the nu- 
clear statistical model, but these theoretical calculations 
are, at present, hampered by large uncertainties in the 
Qf-l-nucleus potential in the astrophysically relevant en- 
ergy range. 

We recently [3] have shown that (n,a) cross section 
measurements on intermediate-to-heavy mass nuclides 
offer perhaps the best means for constraining the many 
parameters defining a realistic a-|-nucleus potential. Al- 
though this approach is showing great promise, one draw- 
back is that the comparison between theory and ex- 
periment as cross sections involves sensitivity to other 
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model parameters in addition to the a-|-nucleus poten- 
tial. A more direct comparison of a strength functions 
extracted from the measurements to theoretical expecta- 
tions should avoid these confounding dependencies and 
allow greater insight into possible improvements in the 
model. 

^''^Sm appears to be the best candidate for such 
a study for several reasons. First, the Q-value for 
^^^Sm(7i,a) is fairly large so the cross section, although 
still very small, is expected to be among the largest 
in this mass range. Second, the natural abundance of 
^'^^Sm is fairly large so the necessary isotopically enriched 
sample is affordable. Third, unlike several other poten- 
tial candidates, there is relatively complete information 
on the parameters for the 100 lowest-energy resonances. 
This is important because reliably extracting a widths 
requires reasonably complete knowledge of the neutron 
and 7 widths for the resonances, and comparison to the- 
ory is most meaningful if the spins and parities of the 
resonances are known. For s-wave neutrons incident on 
^■^^Sm (/■^ = ^), J'^ = 3~ and 4~ resonances can be 
formed in the compound nucleus ^^^Sm. 

The a widths for several ^^^Sm(7i,a) resonances for En 
< 700 eV were determined in previous measurements [4- 
8] of this cross section, and these data were compared 
to statistical model predictions. However, these com- 
parisons were hampered by the fact that the a widths 
for most of the resonances in this region could not be 
measured. We have employed an improved detector to 
make a new measurement of the ^^^Sm(n,a) cross sec- 
tion. In addition to extending the range of the measure- 
ments to much higher energies, this new detector resulted 
in a much improved signal-to-noise ratio in the resolved 
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resonance region so that the a widths for most of the res- 
onances in this region could be extracted from the data. 
This much improved a-width information makes possi- 
ble a more meaningful comparison between theory and 
experiment. 



II. EXPERIMENT 

The experiment has been described elsewhere [3, 9], so 
only the salient features will be given here. The mea- 
surements were made at the Oak Ridge Electron Linear 
Accelerator (ORELA) [10 12] white neutron source. The 
ORELA was operated at a repetition rate of 525 Hz, a 
power of 6 - 8 kW and a pulse width of 8 ns. A 0.76-mm 
thick Cd filter was used to eliminate the overlap of slow 
neutrons from previous pulses and a 1.27-cm thick Pb 
filter was used to help reduce overload effects from the 7 
flash at the start of each neutron pulse. Neutron ener- 
gies were measured via time-of-flight. The detector was 
a compensated ionization chamber (CIC) [13]. Although 
a CIC can have poorer pulse-height resolution than, for 
example, a gridded ionization chamber, it reduces over- 
load effects due to the 7 flash at the start of each neutron 
pulse by several orders of magnitude, allowing measure- 
ments to be made with excellent signal-to-noise ratio to 
much higher neutron energies (500 keV in the present 
case). 

The source-to-sample distance was 8.835 m and the 

neutron beam was collimatcd to 10 cm in diameter at the 
sample position. Two samples were placed back-to-back 
in the center of our parallel-plate CIC with the planes of 
the samples perpendicular to the neutron beam. Hence, 
the cross section was measured over nearly the entire 47r 
solid angle. The samples were in the form of Sm2 03 en- 
riched to 95.3% in ^^^Sm and were 5.0 mg/cm^ thick by 
11 cm in diameter. The ^Li(n,a)^H reaction was used to 
measure the energy dependence of the flux and to nor- 
malize the raw counts to absolute cross section. These 
measurements were made during a run when one of the 
^^^Sm samples was replaced by a ^Li sample. The neu- 
tron energy scale also was calibrated during this run using 
dips in the time-of-flight spectrum caused by resonances 
in the Cd filter and Al vacuum windows in the flight 
path as well as the peak due to the 244.5-keV resonance 
from the ^Li(n,a)^H reaction. A ^Li sample in a separate 
parallel-plate CIC was used as a flux monitor. The most 
recent ENDF evaluation [14] for the ^Li(n,a)^H reaction 
was used in calculating the absolute cross sections. The 
data were corrected for the small background due to the 
spontaneous a decay of ^"^^Sm and for the effects of a 
straggling in the samples. This latter correction (14%) 
was calculated using the computer code SRIM [15]. The 
overall normalization uncertainty of approximately 6% is 
dominated by the uncertainty (± 4%) in this correction 
and by uncertainties (± 3%) in the sample sizes. 



III. RESONANCE ANALYSIS AND RESULTS 

The data were fltted with the 7?.-matrix code SAMMY 
[16] to extract the a widths for resonances in the resolved 
region below 700 eV. Almost all observed [17] resonances 
in this energy range have been assigned [18] as J'^ = 3~ 
or 4~ (s-wave). A radius of 8.3 fm was used in all ^''^Sm 
channels. Both the fact that the sample was an oxide 
and the aluminum backing were included in the input 
files for SAMMY so that corrections could be applied for 
attenuation and multiple-scattering effects in the sample 
and its backing. The resonance energies, spins, parities, 
and neutron and 7 widths from the compilation of Ref. 
[17] , which are based mainly on the work of Refs. [18-21] , 
were used as starting values in the analysis. Only a few of 
the energies had to be adjusted, and one (tentative) spin 
assignment (for the 659.4-eV resonance) was changed to 
fit the data. A radiation width equal to the average for 
^■^^Sm resonances (69 meV Ref. [17, 22]) was used for 
resonances without Tj values in Ref. [17]. The rcsidting 
parameters are given in Table I where, for example, 1.273 
(61) is used to denote 1.273±0.061, etc. Representative 
plots of the data and 7?.-matrix fits are shown in Fig. 1. 

TABLE I: ^*'^Sm{n, a) resonance parameters. 
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The accuracy of the extracted a widths depends on the 
accuracy of the J'', r„, and assignments for the res- 
onances because our measurement technique determines 
only the resonance areas, = gjTaTn/^, where gj is 
the statistical factor (gj = (2J + l)/[(2/ + l){2i + 1)], 
where J, /, and i are the spins of the resonance, ^^^Sm, 
and the neutron, respectively). The uncertainties in the 
resonance areas in Table I are the one-standard-deviation 
uncertainties determined in fitting the data. The uncer- 
tainties in the a widths include additional contributions 
(added in quadrature) from the uncertainties in the neu- 
tron and radiation widths. The uncertainties in the neu- 
tron and radiation widths were taken from Ref. [17]. For 
resonances with unknown radiation widths a "factor of 
two" (± 23 meV) uncertainty was assumed, which should 
be a conservative overestimate of the uncertainty because 
the measured radiation widths [18, 21] show very little 
variability. 

For 23 of the 104 resonances in this region, the fitted 
resonance areas had a relative uncertainty greater than 
70%. In these cases we give only upper limits for the 
resonance areas in Table I equal to the fitted values plus 
the one-standard-deviation uncertainties determined by 
SAMMY. The upper limits on the a widths given in Ta- 
ble I for these cases were calculated from the upper limits 
on the resonance areas using the listed spins and neutron 
and radiation widths. The uncertainties in the neutron 
and radiation widths have negligible effects in these cases 
and hence were not taken into account when calculating 
these a- width upper limits. 



A. Average resonance parameters 

For comparison to statistical model calculations, distri- 
butions or averages of resonance parameters are needed. 
Some of these quantities were determined from a reso- 
nance analysis of the total cross section [21]. With the 
resonance spin and 7-width information from Ref. [18] 
and the a-width information from this work, it is possible 
to calculate these quantities for each of the two possible 
spin states for s-wave resonances and to include the a 
channels as well as the neutron channels. Useful quan- 
tities for comparison to statistical model calculations in- 
clude strength functions, S, average level spacings, D, 
average widths, < F >, and the distributions of widths 
(see Sec. VI for a detailed discussion of these quantities). 
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FIG. 1: ^'"'Sm(n,a) cross section data (points with error bars) 
and SAMMY fit (solid curves). Bars below the data show the 
locations of the fitted resonances. The length of each bar is 
proportional to the square root of the resonance area. 

Care must be taken to account for effects due to missed 
resonances. 



1. Neutron and gamma channels 

The s-wave neutron strength function, =< r° > 
/Z?o, where < > is the average s-wavc reduced neu- 
tron width (r„ = r° X v'^) and Dq is the average 
s-wave level spacing, can be determined from the slope 
of a plot of the cumulative reduced neutron width versus 
resonance energy resulting from a resonance analysis of 
total cross section measurements. This well-known tech- 
nique for determining the neutron strength function is 
relatively insensitive to missing resonances because only 
resonances having small neutron widths are expected to 
be missed in total cross section measurements and the 
slope of the cumulative reduced neutron width versus 
resonance energy is little affected by missing resonances 
with small widths. In Ref. [21], this technique was used 
to calculate 10^ 5o = 4.8 ± 0.5. Using the resonance spin 
information from Ref. [18] together with the neutron 
widths from Ref. [21], strength function plots for the two 



FIG. 2: Strength function plots for the a (top) and neutron 
(bottom) channels. Solid and dashed staircase plots represent 
results for 3^ and 4^ resonances, respectively. The measured 
strength functions are equal to the slopes of these staircase 
plots. 



possible s-wave spins were constructed and are shown in 
Fig. 2. Linear fits to the data in Fig. 2 yield strength 
functions of (4.6 ± 1.0) x 10"* and (4.3 ± 0.9) x 10"* for 
J = 3 and 4 resonances, respectively, in agreement with 
the combined s-wave strength function of Ref. [21]. Un- 
certainties in the strength functions were estimated from 
the number of observed resonances as described in Ref. 
[23]. 

In Ref. [21], the well known technique of determining 

the level spacing from the inverse of the slope of the cu- 
mulative number of resonances versus resonance energy 
(in the lower energy region where missed resonances are 
insignificant) was used to calculate an average level spac- 
ing Do = 5.7 ± 0.5 eV for all s-wave resonances. Using 
this level spacing and assuming that the number of res- 
onances for each spin is proportional to 2J + 1, leads to 
A'' = 54 and 69 resonances by 700 eV for J = 3 and 4 
resonances, respectively (implying Dq = 13.0 and 10.1 
eV for J = 3 and 4 resonances, respectively). These level 
spacings, together with the above strength functions for 
J = 3 and 4, yield < r° >= 5.98 and 4.4 meV, for J = 3 
and 4 resonances, respectively. 
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FIG. 3: Integral distributions of total a (top) and reduced 
neutron (bottom) widths for 3^ (left) and 4~ (right) reso- 
nances in ^'''^Sni. The histograms show the measured number 
of resonances having widths greater than a given width ver- 
sus that width. In an attempt to depict the experimental 
uncertainties, histograms for the a widths axe plotted at the 
measured widths plus and minus their respective uncertain- 
ties. The dashed curves show the expected distributions 
for the average widths, number of resonances, and degrees of 
freedom {v = 1 for neutrons) indicated. 



Neutron widths are expected to obey Porter-Thomas 
[24] distributions having average widths and numbers of 
resonances consistent with the above strength function 
and level spacing determinations. This was found to be 
the case in Ref. [21] for all (sum of J = 3 and 4) s-wave 
resonances. The separate neutron width distributions for 
J = 3 and 4 resonances are shown in Fig. 3. We will 
discuss how well they compare to the expected Porter- 
Thomas distributions in Section V. 



2. a channel 

Our measurement technique is sensitive to resonance 
areas, = gjTaVn/T; hence, missed resonances can 
have small as well as large (e.g. when r„ is small) a 
widths. Therefore, determining the a strength functions 
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FIG. 4: Measured a strength functions averaged over 100- 
eV intervals for = 3~ (X's) and 4^ (triangles) resonances. 
Strength functions for — 3~ resonances calculated with the 
standard (PI) and modified (P2) NON-SMOKER parameters 
are shown as solid and long-dashed lines, respectively, whereas 
corresponding calculations for = 4~ resonances are shown 
as short-dashed and dotted lines, respectively. 



from the slopes of the cumulative a widths versus neutron 
energy (top part of Fig. 2), will likely result in values that 
are systematically small. On the other hand, because 
missed resonances should have a range of a widths, the 
average a widths should be less affected by associated 
systematic effects and strength functions calculated from 
the average widths should be more reliable. There is 
the additional complication that the data indicate that 
the average a widths as well as the a strength functions 
show considerable variations as functions of energy. For 
example, as shown in Fig. 2, there are large steps in 
the slopes of the cumulative a widths versus neutron en- 
ergy. For this reason, the a strength functions were cal- 
culated over different energy intervals. As shown in Fig. 
4, the a strength functions calculated from the a widths 
{S =< T > /D) averaged over 100-eV intervals show 
considerable variation. In addition, as shown in Fig. 5, 
the ratio of the a strength functions for the two different 
s-wave spins states changes dramatically near 300 eV. 
Because it is possible that this effect is caused by incor- 
rect spin assignments and because the strength functions 
as well as their ratio are relatively constant in the 0-300 
eV and 300-600 eV intervals, we list the strength func- 
tions calculated for these two ranges as well as for the 
entire 0-600 eV interval in Table II. We did not include 
the 600-700 keV interval because there are only about 
half as many resonances having firm J'^ assignments in 
this energy range as in the six lower energy bins. We 
compare these values to statistical model calculations in 
Section VI. 

Distributions of a widths for the two s-wave spin states 
are shown in Fig. 3. Alpha widths for transitions to in- 
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FIG. 5: Ratios (J'^ = 3~ to 4~ resonances) of a strength 
functions averaged over 100-cV intervals. Measured ratios 
are shown as solid circles. Ratios calculated with the stan- 
dard (PI) and modified (P2) NON-SMOKER parameters are 
shown as solid and long-dashed lines, respectively. 



dividual final states are expected to obey Porter-Thomas 
distributions, but we were unable to resolve individual 
a groups because our samples were too thick. There- 
fore, the a width distributions are convolutions of Porter- 
Thomas distributions. This will be discussed more fully 
in Section V. 



IV. COMPARISON TO PREVIOUS WORK 

There is reasonably good agreement between the a 
widths we determined and previous work (as compiled 
in Ref. [17]) except that we have many more measured 
widths and hence many fewer upper limits. There are 
however at least two important differences between previ- 
ous results and the resonance parameters reported herein. 
First, the neutron and gamma widths used to extract the 
a widths typically were not reported in previous work. 
Because the a widths extracted can vary substantially 
depending on the values of these other parameters, com- 
parisons to previously reported a widths without this 
information are of limited value. Second, originally [25], 
we followed the example of Ref. [7] and assigned the 
stronger peaks in our data to the same resonances they 
used even though, as the energy increased, our energy 
scale indicated that the peaks actually corresponded to 
the next higher-energy resonances. As a result, our res- 
onance energies became increasing larger than those of 
Ref. [17] as the energy increased. For example, with this 
scheme the resonance we observed at 439.5 eV would cor- 
respond to the resonance in Ref. [17] at 435.7 eV rather 
than to the much closer resonance at 440.2 cV. After a 
recheck of our energy calibration, using the present data 



as well as data from "3]\fd(n,a) [26] and 95Mo(n,a) [9] 
measurements run under the same conditions using the 
same apparatus, we decided that it was not possible for 
our energy scale to be wrong in this manner. Therefore, 
we assigned the resonances we observed to the closest 
ones in Ref. [17]. As a result, many of the spin assign- 
ments as well as the extracted a widths (because the 
widths were calculated using the wrong spins and neu- 
tron and radiation widths) in our preliminary report [25] 
are incorrect, and several resonances with large a widths 
at the higher energies are now assigned to J'^= 4~ rather 
than ,r= 3". 

Further evidence that our spin assignments are correct 
can be seen in the pulse-height distributions for the res- 
onances. As illustrated in Fig. 6, because parity must 
be conserved, 4~ resonances in ^"^^Sm are forbidden from 
decaying to the 0+ ground state of ^''^Nd by a emission 
whereas 3~ resonances are not. Therefore, the pulse- 
height spectrum for a particles from 4~ resonances will 
not include the highest-energy group corresponding to 
decay to the ground state of ^^''Nd. Because our sam- 
ples were relatively thick, we could not resolve the var- 
ious a-particle groups. However, as shown in Fig. 7, 
there is a discernible difference between the pulse-height 
spectra for unambiguously assigned J'^= 3~ and 4~ res- 
onances. As expected, the 3'^= 3~ resonance at 83.775 
eV occurs at larger pulse height than the J^= 4~ reso- 
nance at 18.340 eV. This figure also demonstrates that 
the pulse height spectrum corresponding to the resonance 
at 667.0 eV clearly favors a 4~ assignment for this res- 
onance. Similar comparisons show that 4^ assignments 
also are favored for several other resonances with fairly 
large a widths which had been associated with 3~ reso- 
nances in previous work [7], lending confidence that the 
energy calibration, and therefore the spin assignments, 
of the present work are correct. 



V. POSSIBLE NON-STATISTICAL EFFECTS 

Analysis of previous ^''^Sm(n, a) measurements [4, 5, 7] 
have hinted at the possible presence of non-statistical ef- 
fects. The improved precision and accuracy of the current 
measurements together with improved resonance param- 
eter information [18] makes it possible to perform more 
robust tests for non-statistical effects. 

In nuclear statistical theory, partial widths F associ- 
ated with the decay of compound nuclear states are as- 
sumed to be described by distributions with v degrees 
of freedom: 

where G is the gamma function, x = < F > is 

the average width, and the distribution is characterized 
by a dispersion equal to 2 < F >^ /j/. Neutron widths 
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TABLE II: Neutron and a strength functions from experiment and optical models. 
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10*S(, 






Alphas 
10®5„ 






Experiment 


Optical Model 


0-600 eV 


Experiment 
0-300 eV 


300-600 eV 


Optical Model 
PI P2 


3 
4 


4.6 ± 1.0 
4.3 ± 0.9 


9.2 
9.2 


0.116 ±0.025 
0.092 ± 0.020 


0.171 ±0.052 
0.058 ± 0.018 


0.0560 ± 0.018 
0.127 ±0.039 
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FIG. 6: Energy- level diagrams depicting the ^'*^Sm(n,a)^''*Nd 
reaction [Q = 10.127 McV). Excitation energies are given 
in MeV whereas laboratory neutron energies En arc given 
in eV. The energy scales of the ^''''Nd and '^'**Sm parts of 
the figure differ by a factor of 17 million. = 3~ levels 
in '^''^Sm can a decay to the ground as well as excited (with 
substantially reduced penetrability) states of "'''^Nd. J'^ = 
4" levels in "**Sm are parity-forbidden from decaying to the 
ground state of ^''''Nd, but can decay to excited states. 



as well as partial a widths for specific angular momen- 
tum values to individual final states are expected to have 

= 1, and hence obey Porter- Thomas [24] distributions. 
On the other hand, 7 widths show much smaller fluctua- 
tions because the much larger number of decay channels 
in the 7-ray cascade implies a distribution with many 
more degrees of freedom. In our present experiment, we 
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FIG. 7; Pulse-height spectra from the ^^^Sm(?i,Q) reaction for 
the 18.340-, 83.775-, and 667.0-cV resonances. The spectra 
have been normalized to have equal areas. 



were unable to resolve the individual a groups, and hence 
only the total a widths summed over the various possible 
final states were measured. Because the partial a widths 
fluctuate independently, the fluctuations of the total a 
widths, r„ = X^^ac) are expected to be smaller owing 
to random cancellations of the partial widths. There- 
fore, the distributions of total a widths are expected to 
be narrower {v > 1) than for the partial widths. The dis- 
tribution for total a widths are complicated convolutions 
of partial distributions with v = 1 and different average 
values < Vac >■ Using a Monte Carlo method, it has 
been shown [27] that this convolution does not have to 
be calculated but that instead equation 1, after the par- 
tial width Tac and its average < > are replaced by 
their corresponding total values r„ and < r„ >, can be 
used to describe the distribution of total a widths, albeit 
with an effective degrees of freedom. 



(2) 



where Pc is the penetrability for a particles in channel c. 

In the present case, the calculated [5] degrees of freedom 
are = 1.8 and 2.5 for 3~ and 4~ resonances, respec- 
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FIG. 8: Reduced neutron (top), a (middle), and 7 (bottom) 
widths for '^^'^Sm resonances versus neutron energy. The re- 
duced neutron and 7 widths are from previous work as com- 
piled in Ref. [17] whereas the a widths are from the present 
work. Measured widths with J = 3, 4, or unknown spin as- 
signments are shown as circles, Xs, and squares, respectively. 
a widths for which only upper limits were determined are 
plotted as triangles. 



tivcly. Therefore, it is expected that the fluctuations in 
the total a widths will be intermediate between that for 
neutrons and gammas. However, as shown in Fig. 8, 
although the neutron and gamma widths fluctuate as ex- 
pected, the a widths do not appear to follow the expected 
behavior. For example, instead of being intermediate to 
the gamma and neutron distributions, the a widths show 
the largest fluctuations. Also, there appear to be regions 
of energy in which the a widths fluctuate about as ex- 
pected and other regions of much larger widths and/or 
fluctuations. These effects are evident for both = 3~ 
and 4~ resonances. 

Recasting the data into plots of the cumulative widths 
versus resonance energy, as shown in Fig. 2 further il- 
lustrates that whereas the neutron widths behave as ex- 
pected, the a widths show indications of non-statistical 
effects. The large steps in the cumulative distributions 
of the a widths might be indications of states in -"^^^Sm 
with non-statistical properties. 



Converting the data to the integral width distributions 
shown in Fig. 3 allows a more quantitative comparison 
with theory. In the case of neutrons, the widths are ex- 
pected to follow a Porter-Thomas distributions with av- 
erage widths calculated from the measured strength func- 
tions and level spacing. As shown in the bottom part of 
Fig. 3, there is good agreement between the observed and 
expected width distributions for neutrons indicating that 
the neutron widths behave as expected when separated 
according to spin. The flt to the Porter- Thomas distribu- 
tion for J = 3 resonances can be improved by assuming 
additional missed resonances and hence a smaller aver- 
age reduced neutron width for this spin. This indication 
that there may be slightly more J = 3 resonances than 
the simple 2J + 1 assumption we used is in agreement 
with the relative sizes of the level spacings predicted in 
NON-SMOKER (see Sec. VI) for the two s-wave spin 
states. 

Because missed resonances are expected to have ran- 
dom a widths, simple averages over the measured widths 
should yield good estimates of the average widths and it 
should be possible to compare to the expected distri- 
butions without correcting for missed resonances. These 
comparisons are shown in the top of Fig. 3. For the 
theoretical distributions, average a widths were cal- 
culated from simple averages of the widths determined in 
the resonance analysis and the degrees of freedom values 
were taken from Ref. [5]. In an attempt to depict the 
uncertainties in the Fq, values, two histogram curves are 
plotted in each part of Fig. 3 corresponding to the a 
widths plus and minus their uncertainties. Resonances 
for which only upper limits were determined are plotted 
at the fitted values plus and minus their uncertainties as 
determined by SAMMY. For the 44 firm 3^ resonances 
<Ta > =1.50 /ueV, whereas for the 45 firm 4~ resonances 
< Fq, > = 1.24 /xeV. As shown in Fig. 3, the theoreti- 
cal distributions are substantially different from the data. 
Furthermore, the agreement cannot be improved by, for 
example, decreasing the average widths and increasing 
the number of resonances in an attempt to correct for 
possible effects due to missed resonances. On the other 
hand, for J'^ = 3~ resonances it is possible to obtain 
fairly good agreement between the data and the theoret- 
ical distribution if the resonance with the largest a width 
at 184.76 eV is excluded (resulting in < r„ >= 1.05 eV 
and N = 43). For J'^ = 4~ resonances however, rea- 
sonably good agreement can be obtained only if 3-4 res- 
onances with the largest a widths are excluded. These 
are the same resonances that cause the large steps in the 
top part of Fig. 2. 

One final piece of evidence on the unusual nature of 
the extracted a widths is revealed in a comparison of the 
average widths for resonances for the two different spins. 
In the energy range of this work, emitted a particles from 
the ^^^Sm(n, a) reaction are below the Coulomb barrier, 
so penetrability is a steep function of energy. Therefore, 
the average a width for 3~ resonances is expected to be 
5-10 times larger than for 4~ resonances because the lat- 
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ter are forbidden by parity conservation from decaying 
to the 0+ ground state of ^''^Nd (and hence on average 
have less energy) whereas the former are not. However, 
as noted above, the a widths averaged over ah resonances 
with firm spin assignments are almost equal for the two 
spin states. Furthermore, the a strength functions for 
the two spins are expected to remain constant over the 
range of our resonance analysis given the smallness of our 
range of neutron energies compared to the energies of the 
emitted a particles = 10.127 MeV). However, as 

can be seen in Fig. 5, our data reveal a striking disagree- 
ment with the expectations. We find that the ratio of a 
strength functions for the two spin states changes rather 
dramatically from rj 3 to « 0.5 near En = 300 eV. 



VI. COMPARISON TO THEORY 

Even though the data show signs of non-statistical ef- 
fects, it is interesting to compare the measured average 



resonance parameters to those computed from optical 
potentials. Although this is independent of the Hauser- 
Feslibach approach, using potentials similar to those from 
established statistical models allows a direct test of the 
potentials and an examination of the various compo- 
nents contributing to the calculated cross section, and 
also should yield a better understanding of non-statistical 
signatures. 



A. Definitions 

The statistical model (Hauser-Feshbach approach) [28, 

29] assumes the compound nucleus reaction mechanism 
and a nuclear level density sufficiently high to be able 
to average across it. In this approach, the cross section 
between a target nucleus i and projectile j proceeding to 
exit channel e (i.e. the a channel in our case) is given 
by: 



J 



fc2 (27, + 1)(7,. + 



1)^^^^+^^ Ttot(£,j,^) WUE,J,^). 



r 



(3) 



In this equation, kj is the wave number of the projectile, 
/j and Ij are the spins of the target and projectile, respec- 
tively, and the transmission coefficients (TCs) T{E, ,/, tt) 
and width fluctuation corrections (WFC) Wj^dE, J,7t) 
are defined by [29]: 



T{E,J,7r) 



2tt 



D{E,J,Tr) 



i^jAE)) 



(4) 



and 



Wj,eiE,J,7r) = 



{E)T,,j^dE) 



^ZiE) 



(^'fME)) 



(5) 



The decaying compound nucleus is characterized by the 
average level spacing D{E, J, tt) of states with spin J and 
parity tt at the excitation energy E. Note that in the 
laboratory Tj contains only the transition to the ground 
state of the target whereas T^. contains a sum over all 
possible transitions in the exit channel. The WFC cor- 
relate the incoming and outgoing channels and account 
for pre-equilibrium effects (Wj,e < 1) by rearranging the 
flux into different channels. In practice, they usually are 
implemented by applying corrections to the calculated 
widths. For our case it is important to note the dif- 
ference between widths obtained with or without WFC. 
The WFC arc a model effect. Therefore, widths calcu- 
lated without WFC should be used when comparing to 



directly measured widths. These are the averaged widths 
(F) found in the above equations. However, in the cal- 
culation of statistical model cross sections, the WFC are 
important and must be included. For the {n,a) energy 
range studied here, this strongly affects the neutron chan- 
nel, introducing a further uncertainty in the theoretical 
modelling of the cross section. Nevertheless, the WFC 
are generally thought to be well understood. 

The strength functions S extracted from the experi- 
mental data are closely related to the theoretical TCs; 



SiE,J,TT) 



T{E,J,tt) 
27r 



(6) 



and generally include all energetically and spin- 
algebraically allowed transitions from the compound 
state {E, J, tt) to the final states. As seen in Fig. 6, this 
includes the ao, ai, a2, ■ ■ ■ transitions in the a channel for 
the 3^ compound resonances but excludes the ao transi- 
tion for the 4~ resonances. In the calculation, above the 
first 10 excited states in ^^"^Nd the sum over individual 
states is replaced by an integration over a level density 
in -'^^''Nd. However, in the present case most of the con- 
tribution to the total TC (or strength function) arises 
from the lowest lying states included explicitly. Thus, 
any possible error in the level density in the final nucleus 
is strongly suppressed and does not affect the TCs. 

The TC for a given transition with specified quan- 
tum numbers is computed by solving the stationary 
Schrodinger equation with a given optical potential. 
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Thus, the TCs are sensitive only to the optical poten- 
tial. 

It becomes evident from the considerations above that, 
when possible, it is preferable to compare calculated TCs 
to strength functions extracted from experimental data. 
The TCs are primary quantities depending only on the 
optical potential; thus, this approach makes it possible 
to disentangle the various contributions to the theoretical 
uncertainty. Converting strength functions to theoretical 
widths adds the uncertainty in the compound level den- 
sity p{E, J, tt) = 1/D{E, J, tt). Comparing cross sections 
includes all possible errors in the widths of the different 
channels (optical potentials), the level density, and the 
WFC. 



B. Comparison of strength functions 

Strength functions calculated employing the same 
methods used in the statistical model code NON- 
SMOKER [30-32] are compared to the data in Figs. 4 
and 5 and Table 11. These calculations were made to il- 
luminate the various confounding effects that enter into 
calculations of cross sections and reaction rates as well as 
to ascertain whether it is possible to reproduce the data 
using optical model strength fimctions, especially given 
the indications of non-statistical effects noted above. 

Several differently parametrized potentials are avail- 
able in literature, but for simplicity we limit our inves- 
tigation to one basic shape (Saxon- Woods) and two ba- 
sic parameter sets; the one of Ref. [33] (potential PI, 
standard NON-SMOKER settings) and the other from 
Refs. [34 36] (potential P2). Cross sections for the 
"'^^^Sm(n,a)^''^Nd reaction calculated using potential PI 
are a factor of 3.3 larger than the data [3], whereas calcu- 
lations made using potential P2 are in significantly bet- 
ter agreement (a factor of 1.4 higher than the data) with 
these data as well as data from a number of other reac- 
tions. 

Both potentials use standard Saxon- Woods shapes in 
the real and imaginary parts of the radial potential U : 



l + exp(r^^) l + exp(r^^) 

_ , (7) 

The parameters arc given in Table 111. It is interesting to 
note that the imaginary parts of the two potentials are 
the same. Nevertheless, P2 yields a value closer to the 
measured cross section. At first glance, this seems coun- 
terintuitive as it is commonly stated that the imaginary 
part of the optical potential determines the TC. How- 
ever, it is more correct to state that the TC is given by 
the imaginary part of the wave function which in turn de- 
pends on the relative strengths of the real and the imag- 
inary parts of the potential. 

As can be seen, the calculated strength functions in the 
neutron channel are about a factor of two larger than ex- 
periment for both 3~ and 4~ resonances. It is interesting 



TABLE III: Parameter of the basic Saxon- Woods potentials 



Potential 


V 


rr 




W 








MeV 


fm 


fm 


MeV 


fm 


fm 


PI [33] 


185.0 


1.40 


0.52 


25.0 


1.4 


0.52 


P2 [34-36] 


162.3 


1.27 


0.48 


25.0 


1.4 


0.52 



to note that the standard NON-SMOKER [31] level spac- 
ings {Do = 7.4 and 6.5 eV for J = 3 and 4 resonances, 
respectively) are about a factor of two smaller than the 
measured ones; hence, the calculated average neutron 
widths are in fairly good agreement with the measured 
values. This illustrates the importance of making the 
comparison between theory and experiment as strength 
functions, for if instead average widths were compared, 
the confounding influence of the level spacing might lead 
one to conclude that there was better agreement between 
theory and experiment than there actually is. 

For potential PI, in the i?„ = — 300 cV range, the 
calculated 4^ a strength function is about a factor of 
two larger than measured, but the calculated a strength 
function for J'^ = 3~ is in more serious disagreement (a 
factor of about 6 larger) with experiment. In the En = 
300 — 600 eV range, the 4~ a strength function calculated 
with potential PI is in good agreement with the data, 
but for 3~ resonances the a strength function is in even 
more serious disagreement (about a factor of 17) than it 
was in the lower-energy region. Potential P2 is clearly 
better than PI in predicting the a strength function for 
3~ resonances and for 4~ resonances in the £"„ = — 300 
eV region. However, as shown in Fig. 5, both potentials 
strongly overpredict the 3~ to 4~ a strength function 
ratio, potential P2 yielding a ratio only slightly smaller 
than that calculated with potential PI. 

These comparisons as strength functions indicate that 
the optical a+^'^'^Nd potential requires more adjustment 
than would be surmised from comparisons as cross sec- 
tions. However, it should be noted that the cross-section 
comparisons were made at higher energies {En « 10—500 
keV) where contributions from p-wave resonances are ex- 
pected to become more important. Furthermore, the dis- 
agreement between calculated and actual strength func- 
tions might become smaller at higher energies as has been 
found in other reactions. Nevertheless, it is informative 
that both the calculated neutron and a strength func- 
tions are too large - a fact that could not be surmised 
from comparisons as ^^^Sm(n,a)^'*^Nd cross sections. 

As shown in Fig. 5, the ratio of the calculated a 
strength functions for the two s-wave spin states is con- 
stant over the energy range of resonance analysis. This is 
in stark contrast to the measurements which show a steep 
decrease in this ratio above £"„ = 300 eV. Regardless of 
the potential used, it will never be possible to achieve 
such an abrupt change in the ratio within an optical 
model of strength functions. However, it is possible that 
the abrupt change in the 3^ to 4^ a strength function 
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ratio is due to incorrect spin assignments above En = 300 
eV. Therefore, it is interesting to study if calculations can 
reproduce the measured ratio below En = 300 eV where 
the data should be very reliable. Two conclusions can be 
drawn from a more systematic variation of the parame- 
ters of a Saxon- Woods a potential which we attempted. 
Firstly, the absolute values of the strength functions are 
far more sensitive to the potential than is the 3~ to 4~ 
a strength function ratio. Trying to reduce to calculated 
ratio by only the smaller amount needed to reproduce the 
measured ratio in the region below S„ = 300 eV quickly 
leads to strength functions which are orders of magni- 
tude larger than the observed ones (and hence also to an 
inferior description of the cross section). Secondly, it is 
impossible to obtain a 3~ to 4~ a strength function ra- 
tio smaller than unity with any Saxon- Woods potential. 
In consequence, it is impossible to reproduce the small 
ratio observed above En = 300 eV with any such calcu- 
lation. Thus, from this analysis it appears that a 3~ to 
4~ a strength fimction ratio smaller than unity may be 
an additional indication of a non-statistical effect in the 
data. 

Finally, the NON-SMOKER calculation predicts that 
78% of the ^'*''Sm(n,a)^^'*Nd cross section is given by 
transitions from 3~ states, 12% from 4~ states, and 
10% from other states (higher partial waves). Of the 
3~ transitions, 67% directly populate the ground state 
of ^^^Nd, which cannot be reached from 4~ resonances. 
New measurements in which the various a groups are re- 
solved would be very useful for testing these predictions. 
The relative contributions of the various transitions also 
shows why it is more important for cross section predic- 
tions to reproduce the 3~ TCs than the 3~ to 4~ ratio. 
Therefore, potential P2, which has essentially the same 
3~ to 4~ ratio as potential PI but comes much closer to 
reproducing the 3~ TCs, yields a cross section in better 
agreement with the measurements. 



VII. CONCLUSION 

Comparing theory to experiment as a strength func- 
tions rather than as ^^^Sm(n,Q:) cross sections avoids con- 
founding effects due to the neutron and gamma chan- 
nels, level densities, and width fluctuation corrections 
and therefore can reveal more useful information about 
possible improvements to theoretical models by isolating 
effects due to the a+nucleus potential. Furthermore, sep- 
arating the data into the two possible s-wave spin states 
may yield even more information about the a-f-nucleus 
potential, because a particles from 3~ resonances have 
(on average) larger energies than those from 4~ reso- 
nances and hence they sample a different region of the 
a-|-nucleus potential. Therefore, differences between the 
measured and calculated a strength functions for the two 
different s-wave spin states should be useful for future 
improvements in the a+nucleus potential. An improved 
a-|-nucleus potential would be very useful for astrophys- 



ical applications. 

Interestingly, it is clear from the data presented in Fig. 
8 that the extracted a widths exhibit fluctuations dif- 
ferent from the expected behavior and hint at possible 
non-statistical effects. One striking feature is that the a 
widths exhibit peaks and/or regions of large fluctuations 
as a function of neutron energy instead of the expected 
random fluctuations intermediate to that for neutron and 
7 widths. Perhaps these peaks are a manifestation of a 
nuclear structure effect in the ^^^Sm compound nucleus 
that is not observable in decay channels other than the 
a channel. One important difference between the de- 
cay of the ^"^^Sm compound nucleus into the a channel 
compared to decay into the neutron or 7 channels is the 
large Coulomb barrier that the a particles must over- 
come. Hence, it is possible that the Coulomb barrier for 
a decay could act as a lever arm enhancing the signa- 
ture of nuclear structure effects that are too subtle to be 
observed in other channels. For example, a decay could 
be enhanced for compound states with significant defor- 
mation and for J"^ — A" states of significant collectivity 
(because they decay mainly to the = 2+ first excited 
state of "4Nd). 

Further evidence for possible non-statistical behavior 
of the a widths is revealed when the resonances are sep- 
arated according to spin as shown in Figs. 2 through 5. 
The striking disagreement with theoretical expectations 
shown in these figures depends on reliable spin assign- 
ments for the resonances. The spin assignments used 
in these figures rely on the results of Ref. [18], the good 
match between the energy scale for this reference and our 
work, and the crude pulse-height information from our 
experiment for resonances having large a widths. Below 
approximately 300 eV, the spin assignments used herein 
should be very reliable. Therefore, it is very interesting 
that our exploratory calculations show that, even in this 
energy range, it is not possible to reproduce the observed 
a strength functions as well as the 3~ to 4~ a strength 
ratio at the same time using a Saxon- Woods potential. 

The a-widtli distributions shown in Figs. 2 and 3, as 
well as the striking change in the 3~ to 4~ a strength 
ratio near 300 eV shown in Fig. 5 depend on reliable 
spin assignments above 300 eV. For a number of reasons, 
the spin assignments in this region may not be as reli- 
able, so it would be very useful to make new ^^''Sm(n,a) 
measurements with thinner samples to check spin assign- 
ments, especially for those resonances having the largest 
a widths. In such measurements, resonances having vis- 
ible ao groups to the ground state of ^^''Nd unambigu- 
ously could be assigned as J = 3. The a-particle spectra 
for the few resonances below En — 185 eV that have 
been reported [4] are in agreement with the accepted 
[17] spin assignments, except possibly the 58.130-eV res- 
onance (57.9 eV in Ref. [4]), which is assigned J'^ = 3^ 
but has an almost invisible ao group. New neutron cap- 
ture and total cross section measurements on ^^^Sm also 
could be useful. New measurements with higher resolu- 
tion and sensitivity could reduce uncertainties in the a 
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widths by identifying sonic of the missing resonances in 
the relevant energy range and by providing more 7 widths 
as well as more precise neutron widths. Finally, {n,a), 
(n,7), and neutron total cross sections on other nuclides 
in this mass region should be very useful in shining more 
light on this interesting problem. 
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